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Целью работы является повышение эффективности управления 
транспортными средствами, поэтому помимо синтеза модели 
приведен пример анализа динамики современного пассажирско- 
го электровоза ЭП10. В качестве предмета работы выступают 
математические модели, отражающие динамику вертикального 
движения транспортного средства с учетом нелинейного преоб- 
разования углового движения колеса в поступательное движение 
центра масс колеса. Рассматривается контакт колеса и дороги. 
Представлен синтез модели, отражающей динамику вертикаль- 
ных колебаний транспортного средства. От известных моделей 
приведенная в данной статье отличается тем, что в ней учтена 
динамика углового движения. Таким образом, полученная мо- 
дель связана с координатами состояния углового движения ко- 
леса. Результаты анализа представлены фазовыми портретами, 
отражающими динамику вертикального движения центра масс 
колеса, а также динамику вертикальных колебаний кузова. По- 
лученные результаты позволяют сделать следующий вывод. 
Существует связь между динамикой вертикального и углового 
движения в транспортном средстве. Как следствие, скорость 
углового движения влияет на длину площадки контакта колеса с 
дорогой (рельсом). Результаты могут быть применены при син- 
тезе систем управления транспортными средствами, учитываю- 
щих выявленную связанность. 


Ключевые слова: динамика транспортного средства, площадка 
контакта, вертикальные колебания транспортного средства. 
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Введение. В настоящее время довольно остро стоит проблема удешевления железнодорожных и автотранспортных 
перевозок. Существенно сократить издержки в этой сфере позволит проектирование современных систем автоматиче- 
ского управления транспортными средствами, в контуре управления которыми нет человека. Для современных слож- 
ных динамических транспортных систем принципиальное значение приобретают вопросы построения моделей, опи- 
сывающих контактное взаимодействие колеса и дороги в координатах состояния системы управления [1-5]. Для син- 
теза систем управления транспортным средством необходимо решить ряд проблем —Щ в частности, проблему учета 
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влияния динамики вертикальных колебаний центров масс транспортного средства на тяговые и скоростные характе- 
ристики процесса управления. 
Задача синтеза модели, описывающей вертикальные колебания, не нова. Ее решению посвящено значительное коли- 
чество работ — как о железнодорожной технике [6, 7], так и об автомобильной [8-11]. Отметим, что в этих моделях не 
учитывается взаимовлияние координат вращательного движения колеса и координат, описывающих вертикальное 
движение центра масс колеса и кузова транспортного средства. Вместе с тем, такое влияние подтверждено экспери- 
ментально [10, 11]. Это позволяет сделать вывод о необходимости математического описания данного явления. Итак, 
синтез качественной модели, отражающей указанную связь, является актуальной научной и практической задачей, 
чему и посвящен представленный материал. 

Модель динамики вертикальных колебаний. При построении моделей, описывающих сложную динамику 
систем управления, часто возникает противоречие между требованием к точности модели и ее сложности. Однако для 
качественного анализа основных форм колебаний в транспортном средстве достаточно принять двухмассовую модель, 


которая в упрощенном виде представлена на рис. 1. 
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6) 
Рис. 1. Модель динамики транспортного средства: в координатах состояния (а); 
упрощенная двухмассовая модель (6) 


На рис. | ту, 712 — массы колесной пары и кузова соответственно; ст, с› — характеристики жесткости колеса и 
подвески соответственно; ил, и2 —коэффициенты, характеризующие диссипацию в колесе и в подвеске соответствен- 
но. Горизонтальной пунктирной линией отмечены уровни, соответствующие положению статического равновесия 
центра масс кузова и колеса. Начало отсчета поместим на этих уровнях, а ось 02 будет направлена вверх. Если центр 
масс кузова или колеса находится выше положения статического равновесия, то их координаты 21, 22 > 0, если ниже, 
ТО 21, 22 < 0. С точки зрения анализа вертикальных колебаний транспортного средства, возникающих в нормальной 
плоскости, удобно представлять колесо в виде упруго-жесткой связи. Так, к примеру, делают при решении задач мо- 
делирования в автомобильной технике [8—11]. 

Значение деформации подвески с достаточной степенью точности можно определить как: 

А/ =б 


—2. +2, 


стат.2 





т ы 
где би. = 28 __ статическая осадка пружины подвески под действием силы тяжести. 
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Для построения модели, отражающей динамику вертикальных колебаний центров масс колеса и кузова авто- 
мобиля, используем широко известный подход, предполагающий использование уравнения Лагранжа второго рода 
[4]. Однако при этом примем в расчет влияние времени существования контакта. Для этого с учетом направления осей 
2, 22 введем дополнительную функцию, отражающую эффект «всплывания колеса», который проявляется при увели- 
чении скорости вращения. Подобный эффект достаточно хорошо изучен в динамке металлорежущих станков. Напри- 
мер, в [12, 13] этот эффект объясняется гидродинамической составляющей сил, возникающих в зоне контактирующих 
поверхностей. В рассматриваемом случае учитывается существование контакта не только в режиме относительного 
скольжения колеса, но и в режиме качения. Поэтому «всплывание колеса» можно объяснить тем, что контакт колеса и 
дороги формируется не мгновенно, а в течение некоторого времени. Под временем формирования контакта колеса и 
дороги мы понимаем время, которое будет потрачено на деформацию в нормальной плоскости колеса. Обозначим 
время формирования контакта, которое также является и временем существования контакта, величиной Лё.. Не вдава- 
ясь в глубокий анализ процессов, происходящих в зоне контакта при сближении поверхностей колеса и дороги, отме- 
тим, что существует некоторое эффективное значение Л, при котором сцепление колеса и дороги достигает некото- 
рого максимального значения и при котором деформация колеса в нормальной плоскости останавливается. 

С учетом вышеизложенного и при условии, что ДЁ< Лё,, значение деформационных изменений колеса в нор- 
мальной плоскости можно определить: 


А1, = Он —2 — (О ? (1) 


где (АВ —Щ_ некоторая неизвестная функция, учитывающая влияние времени существования контакта на величину 


и (т, +т,)& 
с 


деформационных изменений колеса в нормальной плоскости, а б — статическая деформация (осад- 


стат.1 


ка) колеса под действием веса кузова и колеса. 
Если ДЕ> Д&, то 


Д1 =б 2. (2) 


стат.1 


При ДЕ < ДЕ. будут справедливы два следующих предела: 
Ни / (АА) > 
41—10 


Вт Г(АЭ — 0. 


АА, 
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(3) 


т ДЕ ы 
С достаточнои степенью точности такое изменение ) можно описать при помощи следующеи функции: 
в, 


(0 — о |= ий * (4) 


Подставим (4) в (1) получим: 
Ав 


АЛ г ба а 2 < (5) 


Для дальнейшего анализа уравнения (5) выразим ДЁ и Д&. в координатах состояния. Для этого рассмотрим за- 
висимость времени существования контакта от длины площадки контакта и скорости углового движения колеса. 

Пусть А — радиус колеса; 9 — угол между двумя прямыми, соединяющими центр колеса с двумя противопо- 
ложными краями площадки контакта; АЛ — величина, характеризующая степень деформации колеса, под действием 
нормально приложенной нагрузки (5); Г, — формируемая длина площадки контакта. 

С учетом принятых обозначений можно определить угловую величину сектора колеса, принимающего уча- 


стие в контакте — 6. Для этого найдем /[: 
Е=2. [4 (2®-А1), (6) 


а затем, используя известные тригонометрические соотношения, определим значение 0: 





А (2В- М 
и т ,) (7) 
© у 


Для железнодорожной техники угол 0/2 достаточно мал. С учетом этого, используя известную теорему о пре- 
дельном значении бесконечно малых величин, будет справедливо: 


9=2 тем. (8) 
[©] 
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Время формирования площадки контакта можно связать с ее угловой характеристикой (8) и скоростью вра- 
щения колеса следующим образом: 


ЖЕ (9) 
© 


где © — угловая скорость колеса. 
Отметим, что контактируют и колесо с дорогой (в железнодорожной технике — с рельсом), и дорога с коле- 


- 0 
сом. Время формирования контакта колеса с дорогои будет определяться как ДЕ= — ‚ а время формирования контакта 
[0 


7Й 
дороги с колесом — как Д! = т 


Е 0 

Для режима торможения будет справедливо Д!= р а для остальных режимов (тяги и боксования) — ДА =—. 
© 

Нелинейная зависимость 6 от АД несколько усложняет наши рассуждения, но (8) легко раскладывается в ряд 


Тейлора. Это позволяет сделать вывод о допустимости линейной интерпретации функции /(А1)в виде: 

ЛА) = К, , (10) 
где коэффициент К» будет определяться углом наклона характеристики / (АЛ) в зоне возможных изменений ДИ! и ра- 
диусом колеса (для К = 0,5 Кь= 14): 

я 2К АЛ 
Ко 
Если известно время формирования контакта Д&. (а мы считаем, что это эмпирически определимая величина), 


4: (0 


то также известно и значение скоростей сх и Г», то есть это заранее заданные константы. 
Преобразуем выражение (5) к виду 


о 


41 =5 ед. (12) 


© 
стат 
Как ВИДНО ИЗ (12), от скорости вращения колеса зависит степень его деформации под действием нормально 
приложенной нагрузки. Это позволяет говорить о взаимосвязанности динамики углового движения колеса (колесной 
пары в железнодорожной технике) и динамики вертикального колебательного движения в нормальной плоскости цен- 


тров масс колеса и кузова. 
С учетом вышеизложенного и направления осей 21, 22 ВЕЛИЧИНЫ деформации колеса и подвески: 


о 


[03 
—2 


АЛ — бе 
А =5 


Уравнение Лагранжа, описывающее динамику взаимосвязанного движения центров масс кузова и колеса в 


12 (13) 


—25 +2. 
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нормальной плоскости, примет следующий вид: 


42, 42) 42 
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Система (14) описывает динамику вертикального движения центров масс кузова и колеса. Отметим, что, если 
угловое движение отсутствует, то последние два члена второго уравнения системы будут равны нулю, и система (14) 
преобразуется к линейному виду. 

Рассматриваемую систему от предложенных ранее [8—11] отличает связь координат состояния, описывающих: 

— вертикальное движение центра масс колеса, 

— динамику углового движения. 

Кроме того, с учетом связи через проскальзывание в данном случае принимается во внимание и поступатель- 
ное движение. 

Результаты численного моделирования. Проанализируем динамику вертикальных колебаний, описываемых 
системой (14). Для этого в качестве параметров уравнений, входящих в систему, подберем величины, соответствую- 
щие современному пассажирскому электровозу ЭП1О [1]. В примере, представленном в этой работе, т> = 50000 кг; 
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с2 = 1.107 Н/м; ц2 = 40000 Н/м/с); ти = 300 кг; с! = 1,6:107 Н/м; ви = 4000 Н/(м/с). Тогда значение б 
б 


= 0,0314 м; 


стат.1 
оитз = 0505 М; Кстав. примем равным 0,99; ох = 40 рад/с (72 км/ч) при Ё = 0,5. 

Для проведения эксперимента по оценке динамики транспортного средства в качестве возмущающего воздей- 
ствия использовались изменения начальных условий по координате 21. То есть моделировалась неровность дороги 
(рельса), при этом скорость вращения колеса совершала колебания относительно некоторого стационарного значения, 
что может соответствовать порывам встречного ветра. Результаты моделирования вертикальной реакции колеса пред- 
ставлены рис. 2, 3. 


42 
4 





21 
Рис. 2. Фазовый портрет колебательного движения центра масс колеса при Г = 72 км/ч 


Реакция на возникающие неровности, представленная фазовым портретом на рис. 2, носит колебательный ха- 
рактер. Взяты значения встречающихся неровностей, не превышающие степень статической деформации колеса 
(5 
верхности возможна потеря контакта колеса с дорогой. Из опыта эксплуатации транспортных систем известно, что на 
неровной дороге транспортное средство должно уменьшать свою скорость. В противном случае возможна потеря 
управляемости. Для проверки адекватности модели мы увеличили скорость. Результаты моделирования системы 
представлены на рис. 3. 


= 0,0314 м). Однако, как видно из рис. 2, при увеличении скорости транспортного средства на неровной по- 
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Рис. 3. Фазовый портрет колебательного движения центра масс колеса при Г = 135 км/ч 


Как видно из рис. 3, с увеличением скорости движения транспортного средства растет вероятность потери 
управления (контакта колеса с дорогой) на неровной поверхности. Таким образом, адекватность принятой модели 
подтверждается. 
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Кузов транспортного средства тоже реагирует на возникающие неровности. Результатом такой реакции долж- 
ны быть колебания центра масс кузова. Эффект «всплывания» колеса при увеличении скорости движения определяет 
сложный характер этих колебаний. С точки зрения проверки адекватности модели, правильно давать приращения ко- 
ординат состояния системы относительно стационарного состояния центра масс колеса (21) и контролировать реакцию 
на такие изменения центра масс кузова (22). Иными словами, речь идет о моделировании зависимости колебаний кузо- 
ва от величины возникающих на пути неровностей. Результаты моделирования вертикальных колебаний кузова 
транспортного средства представлены на рис. 4. 
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Рис. 4. Фазовый портрет колебаний кузова 


Как видно из рис. 4, кузов транспортного средства совершает более сложное колебательное движение, чем 
колесо. Однако в отличие от динамики вертикальных колебаний центра масс колеса, центр масс кузова меньше реаги- 
рует на колебательные изменения скорости движения транспортного средства, то есть колебания затухают. 

Заключение. Представленная в статье модель позволяет численно определять степень влияния скорости 
транспортного средства на динамику вертикальных колебаний центров масс кузова и колеса. Это позволит проектиро- 
вать сложные системы управления движением транспорта. Результаты численного моделирования позволяют сделать 
вывод об адекватности принятой модели. Колебательный характер нарастания скорости транспортного средства обу- 
словливает колебательный характер движения центра масс колеса вблизи равновесного состояния. При увеличении 
скорости движения увеличивается «всплывание» колеса и отмечается некоторый подъем кузова. При синтезе модели 
была допущена возможность полного «всплытия» колеса — как в системе (3). Однако проведенные эксперименты (см. 
рис. 2—4) это не подтвердили, что связано с влиянием подвески. Таким образом, «всплывание» в модели ограничено 
не только весом, но и упругой силой, связанной с деформацией пружины подвески. Отметим, что рассматриваемый 
подход к моделированию динамики связанного движения транспортных систем уже известен [14-17], но предлагаемая 
В работе модель приводиться впервые. 

В качестве дальнейшего развития предлагаемого подхода авторы рассматривают синтез связанной модели ко- 
лебательного движения всего транспортного средства с учетом взаимодействия углового и поступательного движения 
колеса при проскальзывании. 
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